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Abstract. 

We calculate the variation with temperature of the vortex free energy in D=2+l SU(2) 
lattice gauge theories. We do so both above and below the deconfining transition at 
T = Tc. We find that this quantity is zero at all T for large enough volumes. For 
T < Tc this observation is consistent with the fact that the phase is linearly confining; 
while for T > Tc it is consistent with the conventional expectation of 'spatial' linear 
confinement. In small spatial volumes this quantity is shown to be non-zero. The way 
it decreases to zero with increasing volume is shown to be controlled by the (spatial) 
string tension and it has the functional form one would expect if the vortices being 
studied were responsible for the confinement at low T, and for the 'spatial' confinement 
at large T. We also discuss in detail some of the direct numerical evidence for a non- 
zero spatial string tension at high T, and we show that the observed linearity of the 
(spatial) potential extends over distances that are large compared to typical high-T 
length scales. 



1 Introduction 



The idea that vortices are the relevant degrees of freedom for understanding con- 
finement in non-Abehan gauge theories is an old one Consider, for example, an 
arbitrary SU(2) gauge field, {^^(x)}, in three dimensional Euclidean space-time. Now 
suppose that we perform the following gauge transformation {A} {A^}- Choose a 
space-like plane and choose a point Xc in that plane. Consider any closed loop that 
encloses that point. We perform a gauge transformation such that, as we go around 
that closed loop, the gauge transformation varies from 1 to a non-trivial element of 
the centre of the group: in this case z = —1. Clearly this gauge transformation is not 
single-valued; nonetheless the gauge potential is, because it is in the adjoint represen- 
tation. However any field in the fundamental representation would acquire a factor 
of —1; in particular, the Wilson loop defined on such a closed curve, c, will acquire a 
factor of —1 relative to its value prior to the gauge transformation 

Wc[A<^] = Tr el^"-^' = -Wc[A] (1) 

We extend this construction to all space-like planes, so that the point Xc traces a 
line and we shall refer to this as the vortex (world) line. We can easily generalise 
the construction so that this line may be closed rather than infinite. If the closed 
curve along which the Wilson loop is defined links once with the vortex line, then the 
Wilson loop will acquire a factor of —1. If it links times then it acquires a factor 
of (—1)^. One can consider the field {A^} to be the original field {A} with an added 
ultraviolet centre vortex, whose core coincides with the vortex line. Along this line the 
gauge transformation is singular and the gauge potentials are ill-defined. Since such a 
vortex corresponds to a gauge transformation everywhere except along a line, we can 
readily define fields with any number of such vortices (since the problematic vortex 
intersections will have measure zero). And it is then straightforward to show, as we 
later shall, that if the vacuum contains an uncorrelated gas of arbitrarily long vortices 
then the average value of a Wilson loop bounded by a planar curve c, will decay as the 
exponential of its (minimal) area, A[c\, 

(I^e) ocexp{-a^[c]} (2) 

with the resulting string tension, a being proportional to the average number of vortex 
lines threading each unit of area. 

When the theory is regularised, such an ultraviolet vortex will become well-defined 
but will have a divergent action density along the actual vortex line. For example, 
consider a cubic lattice with a plaquette action. If we place the vortex line so that it 
avoids all sites and links, then the field is clearly well-defined everywhere. However 
the plaquettes threaded by this line will have their signs flipped, and this will corre- 
spond to an action density along the vortex line, that diverges as we approach the 
continuum limit. Quantum fluctuations and renormalisation will presumably lead to 
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the line spreading into a tube of finite width, with a non-zero but finite action density 
within this vortex 'core'. Indeed, at the classical level we would expect this core to 
spread to fill the whole box: in which case such a vortex ceases to be a useful degree of 
freedom. However it may be that quantum fluctuations induce an effective potential 
for these vortices such that the core is stabilised with a width given by some character- 
istic length scale. Moreover if this potential is such that these vortices also condense 
into the vacuum, then they will produce linear confinement between particles in the 
fundamental representation. But all this is, of course, not guaranteed; if, for example, 
the vortices are themselves confined, then they will not disorder large Wilson loops. 

Whether such vortices exist and, if so, whether they condense, are thus the important 
dynamical questions for this approach to confinement. They are also difficult non- 
perturbative questions that have begun to be addressed, during the last two years, in a 
number of lattice Monte Carlo studies. (See, for example, ^ and references therein.) 

Linear confinement corresponds to an area decay, as in (^, for arbitrarily large 
Wilson loops. If a vortex is to disorder such a loop it must pierce the minimal area an 
odd number of times and this implies that the vortex line must be correspondingly long. 
(Unless it is close to the perimeter.) That is to say, one needs a condensate of arbitrarily 
long vortex lines. Actually this is only so at T = 0. At high T the requirement becomes 
far less severe. A vortex that closes upon itself through the periodic temporal direction, 
is only of length 0{1/T) and yet can disorder spatial Wilson loops that are arbitrarily 
large. (It will not disorder Wilson loops that lie in a space-time plane, but that is to be 
expected in the high-T deconfining phase.) If vortices exist at all then there will surely 
exist a finite density of such periodic vortices in the high-T (Euclidean) vacuum. That 
such short periodic vortices might be the dynamics underlying the apparent area decay 
of spatial Wilson loops at high T is an old idea. However if the interactions between 
these vortices were such that they totally screened each other, or were confined, then 
they would not disorder large spatial loops. That something like this might be the 
case has been claimed by a recent study ^ of the properties of the magnetic centre 
symmetry at low and high T. This study also suggests that the high-T spatial string 
tension is only apparent, and that on appropriately large length scales the effective 
spatial string tension will vanish. 

The present work was originally motivated by these @, |^ unconventional claims. 
We set out to perform lattice Monte Carlo calculations to determine the properties 
of the short periodic vortices at high T and to determine whether the spatial string 
tension was indeed non-zero. The claims made in ^ are strongest for D=2-|-l gauge 
theories and so these are the theories we consider. We work with the simplest SU(2) 
non-Abelian gauge theory. Our work has, so far, involved a careful reanalysis of some 
earlier, partially unpublished calculations of the spatial string tension as a function 
of T, together with some simple probes of vortex properties. The results are striking 
and, if taken at face value, appear to disagree with the suggestion of |Q, ^ concerning 
the behaviour of vortices and spatial confinement at high T. Rather they provide some 
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quantitative evidence that vortices at low T and (timelike) periodic vortices at high 
T have a very similar dynamics. Our work may also be viewed as complementing 
recent calculations concerning the vortex free energy in D=3+l low-T SU(2) gauge 
theories. 



2 SU(2) lattice gauge theory in 2+1 dimensions 

We work on a cubic lattice with gauge fields represented by SU(2) matrices, U^{n), on 
the links of the lattice. We use the standard plaquette action 

= E(l - \TrU,) (3) 

where Up is the ordered product of link matrices around the plaquette labelled p. This 
action appears in the partition function as 



j\{dUieM-(3S} (4) 



where 

/3 = A- (5) 

This relation holds in the continuum limit and we use it to define at arbitrary 
values of the lattice spacing a. (Note that in 3 dimensions the gauge coupling g"^ has 
dimensions of mass.) It is clear from (^ that the continuum limit a ^ corresponds 
to /3 — *• oo. This refiects the fact that the theory becomes free at short distances, just 
as in 4 dimensions. 

To study the theory at finite temperature we use a. ■ Lt lattice, where Lt is the 
number of sites in the periodic 'temporal' direction, and L ^ oo in the thermodynamic 
limit. The temperature of the system is given by 

aT = 1. (6) 

This lattice theory has been studied in detail in recent years and we summarise 
here some features that will be useful for us later on. 

2.1 the confining phase at low T 

The theory is linearly confining at low T. The numerical evidence for this is no less 
clear-cut than in 3+1 dimensions. The T = string tension ^ can be accurately 
parameterised as 

1 r 0.945 l.lOl 

av/a=-|l.337+-^ + -^| :/?>3.0 (7) 
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with an error of about ±1%. We can use @ to convert (0) into an expression for the 
dimensionless ratio \J(y j in terms of powers of ag^ . When we need a value of ai/a, 
or of ^/cr /g"^, at a value of /? at which direct calculations do not exist, we shall use 

©• 

An alternative physical mass scale for the theory is provided by the mass gap, mo, 
which, here, turns out to be the mass of the lightest J^^ = 0"*"^ state. An accurate 
parameterisation for the ratio of mc to v^a is given by 

!!i2 = 4.718 - ,,3> 3.75 (8) 

with an error of < ±1% for (3 > 6 and somewhat larger at lower (3. Using and (|^) 
we can convert (H) into an expression for amc or mc/g'^- 

2.2 the deconfining transition at T = Tc 

As the temperature of the system is increased, the theory undergoes a second order 
deconfining phase transition at T = T^. (See e.g. ||, H].) The values of in can be 
accurately parameterised by 



which can be translated into a critical value of L 



Lm = Y^iP-0.37). (10) 

We will use (|^) and ( p!OD to express the value of T = 1/aLt in units of Tc, at the value 
of (3 at which we are working. 

2.3 the deconfined phase at high T 

In the deconfined phase, T > a Wilson loop defined in a space-time plane will 
not decay exponentially with its area. However it has long been believed that purely 
spatial Wilson loops do continue to show an area decay as in (^. That is to say, in 
the deconfined phase one has 'spatial' linear confinement with a corresponding non- 
zero string tension, cr{T), that will depend on the temperature. A number of numerical 



studies have served to confirm this expectation (see for a recent example), although 
it is probably fair to say that the effort put into this question has been very much less 
than that put into the study of confinement at T = 0. This perhaps leaves room for 
the suggestion made in Q that there is in fact no spatial confinement at high T and 
the apparent spatial decay is illusory. 

Let us consider how this might be. The straightforward possibility is that (accurate) 
numerical calculations have not been performed on distances that are larger than the 
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relevant dynamical length scale, and so the apparent area decay simply reflects an 
approximate 'accidental' linearity in the spatial potential at small/intermediate dis- 
tances. An 'accidental' linearity is far from impossible. After all, the approximate 
linearity observed in the T = potential extends down to distances well below the 
characteristic distance scale ~ 1/m at which one would have expected a flux tube to 
have started forming. Thus, the usual heavy-quark potential may be regarded as being 
'accidentally' linear at small distances. What are 'small distances' here? If T is close 
to Tc then 

^ (11) 



might be a relevant distance scale. At higher T (see |TT[| for a review) we expect 
perturbation theory (in the dimensionless coupling g'^/T <^ 1) to be good, so that 
dynamically generated masses are ~ g^T, and hence a typical dynamical length 
scale will be 

= i/gT'/'. (12) 

However this ignores the infrared divergences in perturbation theory. For example, 
the Debye screening mass m/j is infinite at one loop. Fortunately direct numerical 
calculations of m^) exist |TTI] and these provide a reference scale 



= l/rriD. (13) 

We will consider these scales when we come to discuss numerical results on the spatial 
string tension. 

An alternative and more subtle possibility is provided by the actual situation in 
QCD and similar theories where non-Abelian gauge fields are coupled to fundamental 
matter fields [|1^. Here direct numerical calculations of the potential typically find 
it to be linearly rising even beyond the distance at which one expects the string to 



break, and it requires a sophisticated mixing analysis to 'see' the string breaking 
In such theories the breaking of the string is analogous to the decay of a resonance. 
The string state exists, but as an unstable state rather than as an asymptotic state. 
In QCD the decay width is small and so the overlap of a good string operator onto 
the decay products will be too small for the breaking to be visible in a calculation of 
limited accuracy. This scenario requires an underlying confinement mechanism that is 
weakly broken, as in QCD. If there is no confinement to start with, as in weakly coupled 
D=3-|-l QED, then this scenario is irrelevant and there is no sign of any linearity in the 
potential. While we are not aware of any underlying high-T spatial confinement that 
would make the QCD example relevant to the discussion of 0, |^, if it were to be so, 
then the fact that the string breaking is a weak effect should make it easy to identify 
states into which it decays. For example, in the case of QCD it is into a light qq pair 



plus the static sources. Then a mixing analysis [|I2[ could be performed to demonstrate 
the string breaking. 
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3 The spatial string tension for T > Tc 



In this section we shall analyse one particular calculation of the spatial string tension 
in the deconfined phase. We shall look closely at the distance scales over which we see 
linear confinement and ask whether these distances are 'large' or 'small'. We shall also 
discuss the values and behaviour of cr(T), since, later on in the paper, we shall wish to 
compare our results with theoretical expressions that incorporate such a spatial string 
tension. 

We shall need to discuss the results of § in some detail for two reasons. The first is 
that some of what we present here has not been previously published. The second is 
that P] was not overtly interested in finite temperature physics and a careful translation 
of notation is needed. The relevant calculations in ^ were performed on Li x L2 x L3 
lattices. The Polyakov loops were along the 2 direction and hence the periodic fiux 
tubes onto which they projected were of length L2. The correlations were taken in the 
1 direction. The value of L3 was varied and fiux tube masses were obtained for various 
values of L3 and L2. In L3 was labelled L± and was referred to as a transverse 
spatial direction: the point of the calculation being to estimate the intrinsic width 
of the confining fiux tube. If we choose instead to call Lt = L3 the Euclidean time 
direction then what is calculated in is the mass, Mp, of a spatial confining fiux tube 
of length L = L2 at a temperature of aT = 1/Lt. From this we can extract cr(T) in 
the standard fashion: 

aMp = a^a{T)L - ^ ■ y (14) 

where the second term on the RHS is the universal string correction for Polyakov loops. 
These calculations used Polyakov loops separated in the 1-direction and individually 
smeared in a 2-3 plane. The smearing, however, is only a device for improving the 
efficiency of the calculation. The lightest fiux tube mass is the same as one gets with 
unsmeared Polyakov loops. And this should give the same string tension as one gets 
with a large Wilson loop in the 1-2 plane, once the different string corrections have 
been taken into account. 

3.1 linear confinement? 

The calculations in were performed at /5 = 4.5 and aX [3 = 9.0. What evidence do 
we have for linear confinement? Linear confinement will, of course, be refiected in an 
approximate linear rise in the fiux tube mass, Mp, with its length, L. In Fig.|l]we show 
how, at /? = 4.5, this mass varies as the fiux tube length varies from L = 3 to L = 12. 
We show this for two different values of Lj. The first is for Lf = 2 which corresponds 
to a temperature of T = l/2a = 1.33 x Tc (using (H) and (p!0D), which is well in the 
high-T deconfined phase. The second is for Lf = Q which is well in the low-T confined 
phase. We see an approximately linear rise throughout this range of L, for both values 
of T. 
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Consider, in more detail, the Lt = 2 high-T calculation. The accurately calculated 
values are for L < 8. (The L = 12 mass is so heavy that the errors become large.) Is a 
length of 8a, at /5 = 4.5 and at T = 1.33Tc, large or small? Now the dynamical scales 
defined in ([Tl|), (p^) and ( [T3| ) are, here, — 2.7a, — l-5a and ~ 2.7a. (For the 
last relation we have used the value arrin = 0.373(5) which was calculated in p for 
P = 4.5 and aT = 1/2, but which was not published.) Compared to these length scales 
we can confidently claim to have seen linear spatial confinement over 'large' distances. 

Of course the P = 4.5 calculations correspond to a rather coarse lattice spacing. At 
/? = 9, on the other hand, the lattice spacing is well into the weak coupling scaling 
region and here there exist calculations |^ of flux tube masses for two lengths, L = 8 
and L = 16. One again observes an approximate linear rise at both low and high T. 
For example, at = 4 (where T ~ lAT^) one finds aMp{L = 8a) = 0.292(3) and 
aMp{L = 16a) = 0.677(6). If we evaluate the previously discussed high-T dynamical 
length scales at this (3 and Lt, we find — 5.6a, C,t — 3a and C,d ~ 5a. (For the last 
relation we do not have a directly calculated value of rriD and have instead interpolated 
between 'neighbouring' values calculated in ||ll| and unpublished work from 0.) So 
here too we can claim to have seen linear spatial confinement on reasonably 'large' 
distance scales. 

From the above we feel able to say that the spatial linear confinement we observe is 
not an accidental linearity obtained over small distance scales. On the other hand, the 
possibility of a weak string-breaking is something we cannot address without a more 
specific scenario for the string-breaking dynamics. 



3.2 (j(T) 

We take from ||^ the values of the flux tube mass, Mp, at /3 = 9, for various values 
of the temperature aT = l/Lf and for a flux tube length of L = 16a. We then use 



(14) to extract the corresponding string tensions, cr(T), which we plot in Fig.^. We see 
that o"(T) oc T for T > Tc, just as one expects from the simplest dimensional reduction 
arguments. (This observation, in a different notation, was already made in and a 
simple dynamical reason for it was given there.) More specifically we find: 

cr(T) ~ OAOg^T T^T : T > T^. (15) 

At = 4.5 only the Lt = 2 calculation is in the deconfined phase, so we cannot perform 
a similar calculation there. However if we assume a linear rise in T then we obtain a 
coefficient of ~ 0.44 rather than the 0.40 in ([l5|). The size of the difference is consistent 
with being a 0(1/(3) lattice correction and, if so, this implies that the continuum value 
of the coefficient will be ~ 0.36. Clearly, a more detailed calculation to control the 
continuum limit would be useful. 

In a later section we shall need a particular value of the string tension: 

aV(T = l/4a) = 0.0443(4) : /? = 9.0. (16) 
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4 Vortices and boundary conditions 



4.1 centre vortices on an infinite lattice 

Consider an arbitrary gauge field, {U^{n)}, on an infinite lattice. Let the sites be 
labelled by integers {nx,ny,nt). We can add an ultraviolet vortex to this field by the 
following transformation on the fields, U —> U' . Choose some particular integers n^, 
which will label the spatial plaquette threaded by the vortex line. We take C// = Ui for 
all links / except that 

Uy{nx, Uy, nt) = -Uy{nx, n^, rit) for rix < nl and Vrit. (17) 

It is easy to see that this transformation fiips the sign of the spatial (x, y) plaquettes 
emanating from the sites {nl,ny,nt), for all values of rit, but that it leaves all other 
plaquettes unchanged. (Since —1 is an element of the centre of the group and so 
commutes with all Ui.) Moreover if we consider an arbitrary spatial Wilson loop, it is 
either unchanged or its sign is flipped. The latter happens if and only if the Wilson loop 
encircles a plaquette whose sign has been flipped. This has all the expected properties 
of a vortex. It is an ultraviolet vortex because the vortex core, where the action has 
changed, is a string of elementary plaquettes. The change of sign of these plaquettes 
means that their average action has increased; and that the action density along the 
string will diverge in the continuum limit. The associated Z2 gauge transformation has 
also been localised onto an ultraviolet length scale: a single link. Such a vortex may be 
suppressed by its action, but the measure of the transformed field is clearly the same 
as that of the original field. 

In the above construction there is a half-plane of links that acquires a minus sign. 
This half-plane can be smoothly deformed while keeping the edges fixed. On an infinite 
lattice one edge is invisibly far, and we can deform a set of = y negative links heading 
off to n-r = — 00 to, say, a, set oi /i — x negative links heading off to = —00. On a 
finite lattice this will no longer be possible. 

It is also clear that we can introduce any number of these ultraviolet vortices, with 
vortex lines at any positions. By allowing n^,n^ to vary smoothly with rit we can 
introduce vortices whose world lines are not straight. A suitable pair of such vortices 
can produce a vortex line that is a closed loop in space-time. Thus an ensemble of 
ultraviolet vortices may consist of closed loops of arbitrary sizes and shapes. 

When quantum fiuctuations are included, we expect that if vortices do exist then 
the core will acquire a size O(^) where is an appropriate dynamical length scale. 
(And that in any case the gauge transformation will typically be spread over distances 
greater than one link.) However any such field can be reached by a smooth deformation 
of the ultraviolet vortex we have defined above, and so, for the 'kinematic' part of our 
discussion, it suffices to focus on the latter. 
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4.2 twisted and periodic boundary conditions 

From now on we shall be in a finite spatial volume. Suppose that we are on a torus 
so that the gauge potentials are periodic in x and y, with periods L^. and Ly in lattice 
units. If a gauge field {Ui} is periodic then it is obvious that when we add an ultraviolet 
vortex, as in ([T7|), then the transformed field {U[} is not periodic. However we can 
add two such vortices, located at say n^^, n^^ and n^^, n^^, since we can deform the two 
strings of flipped links so that they join the two vortex cores and there is no string of 
flipped links heading off to n^. = — oo. A Wilson loop that encircles both vortices will 
not flip sign but one that encircles only one, will. Thus a periodic gauge field possesses 
an even number of such vortices. 

To obtain fields with an odd number of vortices we need to impose twisted boundary 
conditions. We shall do this in the usual way, keeping the boundary conditions periodic 
but transforming the plaquette action in (|^) as follows (see |TT| for a review): 

St^ = E (1 - ^TrU,) + E(l + IttU,,) (18) 
p¥=p' p' 

where the set of plaquettes p' is specified to be the set of (x, y) plaquettes emanating 
from some specific spatial coordinates n^^,n'^ at all values of nj. That is to say we 
flip the sign of a particular string of spatial plaquettes in the action, the string closing 
on itself through the periodic temporal boundary. That this is equivalent to the usual 
plaquette action with twisted boundary conditions can be shown by making a singular 
centre gauge transformation followed by a suitable redefinition of field variables. 

Consider a field {U} that is typical of those that one obtains with the usual plaquette 
action (^. Now consider the same field, but with respect to the twisted action. The 
action density is unchanged, except along the line of twisted plaquettes where its sign 
is flipped. This is just what we would expect if we had added an (ultraviolet) vortex 
line so as to thread this line of plaquettes. That is to say, we appear to have added a 
single vortex to the field {U}, something that we were previously unable to do. 

This interpretation becomes clearer if the vortex is not superimposed upon the 
twisted plaquettes, but is added to {U} so that it threads some other line of pla- 
quettes, labelled by n^,n^. We do this, as usual, by multiplying by —1 a string of 
links running off from n^,n^. Now, instead of running this string to the boundary 
(which would lead to problems with the periodic boundary conditions) let us run it 
up to n^^,n^^ where the twisted plaquette lies. This will flip the value of the twisted 
plaquettes. So now the action density of this field with a twist is the same as that of 
the original field U without a twist, even at the location of the twist, except that the 
plaquettes at the location of the vortex core have been flipped. This is exactly what 
we got when we added a single vortex to a field in an infinite volume. Moreover, any 
Wilson loop enclosing the vortex, but not the twist, will acquire a factor of —1 relative 
to its value in the field U. We may think of the twist as adding a vortex source to the 
system. Alternatively we note that if we choose to trade the twisted action for twisted 
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boundary conditions, then we obtain a string of flipped links from the vortex to the 
twist and then on to the boundary. But now this presents no problem because if the 
field {U} was periodic, the field with a vortex added satisfies the twisted boundary 
conditions. Thus the ensemble of fields obtained with a twisted action is just like the 
ensemble from the simple plaquette action, but with a vortex added or subtracted. If 
periodic field have even vorticity then twisted fields have odd vorticity. 



4.3 vortex free energy, Fv 

Consider a system at temperature T = 1/aLt. The partition function and free energy 
are related by 

exp|-^l = Z. (19) 



T, 

Let Ftw be the free energy of the system with a twist, and Fnt of the identical sytem 
with no twist. We define the vortex free energy to be the difference 

Fv = Ftw — Fnt- (20) 

Note that this will not be the free energy of a single vortex, unless the usual vacuum 
contains no vortices. It is the free energy difference between a system with an odd 
number of vortices, and a system with an even number. 
From (pop, ([T9| ) and (^ we readily obtain 



-( 










df3^ 






-( 





{Stu.) - {Snt) (21) 

where Stw is defined in (|18|), Snt is the usual untwisted action defined in (^), and we 
transform the derivative with respect to (3 into one with respect to T using the fact 
that the former is performed at fixed Lj. 

If F„ 7^ in the thermodynamic limit [L^^Ly oo) then the vortices certainly do 
not condense in the vacuum and will not contribute to a linearly confining 'spatial 



potential'. This is known to occur for vortices in a space-time plane at high T 
In this paper we calculate not F^ itself, but rather d/dT{Fy/T) using (^Tj). We shall 
assume that if the latter is zero, then so is the former (since we see no reason for Fy to 
be exactly linear in T). A direct calculation of F„ is possible but less straightforward 



4.4 action of a single vortex 

A single ultraviolet vortex imposed on a gauge field with f// = 1 V/ will clearly have an 
action 

5-, = 2Lt = 4^. (22) 
al 
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If the vortex is fat, so that the vr units of flux are smoothly spread over Ay plaquettes, 
then 

^. = 5:(l-cose,)^^^L,. (23) 

If we are working with fluctuating gauge fields at some value of (3 then in general we 
expect the vortex action to be modified by a lattice renormalisation factor: 

SM = Syioo)-Zy{P). (24) 

We note that if, for example. Ay ~ (^d/o)^ = l/{a? ■g'^T) then [3Sy is independent of 
a or T. On the other hand, if the vortex size is not stabilised by quantum fluctuations, 
it will spread over the whole spatial volume, L^, so that 

5.^y-§-^.(/3)'^0. (25) 

In this case the vortex merges into the Gaussian spin-wave contribution and ceases to 
be a useful degree of freedom. This is indeed what occurs in U(l) lattice gauge theories: 
the flux is conserved and hence comes in closed loops, but these do not come in vortices 
of finite area. Although Wilson loops are still disordered, the effect is Coulombic rather 
than linearly confining. 



4.5 domain walls and spacelike vortices at high T 

To give substance to the above discussion of vortices, it would be useful to identify 
a situation where the vortex becomes a real and easily visible excitation, just like a 
magnetic monopole in the Higgs phase of the Georgi-Glashow model. It turns out that 
such a situation exists |Tl[]. Consider a spatial volume L x Lj_ at a temperature T 



close to zero. Suppose also that L is large. If one decreases Lx_ from large values, 
one finds that there is a phase transition at some = L^(/3). If in addition the 
system possesses a spatial twist, then, for L±_ <^ L\{(3), one finds that there is a single 
vortex whose world line extends through time, and which close upon itself through the 
periodic temporal boundary conditions. Moreover the core of this vortex extends right 
across the short direction. In the longer L direction it has a size oc y/{l3L^). As 

3 

L±_ I the vortex action density grows a.s Sy oc 1/Lj_ and the probability of any virtual 
vortex quantum fluctuation rapidly vanishes. Thus, the special spatial geometry drives 
the system into an effective vortex Higgs phase, in which the vortex appears 'squeezed' 
by the short spatial direction. Apart from this distortion, it is precisely a vortex of the 
kind we have been discussing. 

In [|lll the context of these objects was (apparently) quite different, and we need to 
provide a brief translation. Consider the system described in the previous paragraph. 
Rename the time direction to be a spatial one and the short L± direction to be the 
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temporal one. Then the phase transition of the previous paragraph is simply the 
deconfining transition, Tc = l/{aL']_), and the squashed vortex is the 'domain wall' 
separating two high-T Z2 vacua. (Or, more precisely, it is the world sheet of the string- 
like boundary that separates the two Z2 vacua.) The Polyakov loops on either side of 
the vortex have a relative minus sign: just as one would expect for a centre vortex. 
It was briefly pointed out in [|11[| that these domain walls have a simple interpretation 
as real vortices in a space-time plane. Renaming axes it becomes a spatial vortex as 
discussed in the previous paragraph. 



5 Vortices and confinement 

In this section we shall assume that centre vortices exist at high T and that they are 
the fluctuations which produce the spatial string tension. We assume (reasonably) that 
the relevant vortices are the short ones whose world lines close through the temporal 
boundary and whose length is therefore 1/T. These vortices will have a typical area 
Ay, and a corresponding action, 5'^,, as in (|23|J2^ ). Note that we are talking here of the 
vortices in the Euclidean field configurations, rather than the vortices that may appear 
in the eigenstates of the Hamiltonian. Thus may depend on T. 

In the following calculations we shall assume that the vortices are dilute, since it 
makes the argument simpler. Whether this is actually so is something we shall return 
to in the next Section when we discuss our numerical results. As long as vortices are 
uncorrelated at sufficiently large distances, our conclusions will not be qualitatively 
altered by vortex correlations which, to a first approximation, can be absorbed into a 
renormalisation of the effective vortex action. 

We shall first show how large Wilson loops acquire an area decay in the presence of 
such a gas of vortices. We shall ask what the observed T dependence of <j(T) teaches 
us about the properties of these vortices. We shall then show how the average action 
with twist approaches the periodic action as the volume— >■ 00, and how this approach 
is directly related to the string tension generated by the gas of vortices. This parallels 
the well-known results obtained in [|1| for T = 0. 



5.1 vortices and the spatial string tension 

Consider a spatial Wilson loop defined along a curve c that bounds a minimal area A 
(in lattice units). Let us calculate its value within a simple gas of vortices. Suppose 
that the area A is pierced by vortices. Assume they are uncorrelated so that the 
probability is given by a simple Poisson formula 

P{n,;A)= -e-^"-^ (26) 
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where /xq is the probabihty of a vortex piercing a unit area. Now each vortex contributes 
a factor —1 to the value of the Wilson loop, so in this simple picture 



(W,) = {Tr el""-'') = ^(-ir-PK;^) 

nv even ri„ odd 

= e-2^°-^, (27) 



where the last relation follows from using (|26|). That is to say, the gas of vortices 



0- 



produces linear confinement with a string tension a a = 2/i 

Since we are at high T we expect the relevant vortex world lines to be closed through 
the short temporal boundary. If their cross-sectional area is (in lattice units) then 
we expect /iq in (^) to contain a factor ~ as well as a factor containing the action 



aV(r)=2/io ~ Ae-/55" (28) 

~ Ae-l5;7-^-^-(^). (29) 
Ay 

One uncertainty here is in the integration over spatial translations of the vortex. This 
will produce, in addition to the factor l/A^, a factor, presumably 0(1), that depends 
on the details of the collective co-ordinate calculation. There will also be a factor from 
the fluctuations around a vortex. As we remarked earlier, a natural length scale at 
high T is provided by = 1/y/g'^T. If we assume that Ay oc (^d/*^)^ then we find 
from (^) and ( P7| ) that the string tension is 

aV(T) = 2/io = ca^g^T (30) 

where c may have a weak (3 dependence from the Zy{f3) factor. Thus we naturally obtain 
a spatial string tension that increases linearly with the temperature (as one would 
expect from naive dimensional reduction). Conversely the observation that o"(T) oc T 
implies that Ay oc 1/T and hence that the weighting exp{—pSy) is independent of T. 



5.2 AS and vortices 

What do we expect for the action difference 

AS = (Sty,) - {S^t) (31) 

that appears in (|2T|)? The only difference between having twisted and periodic bound- 
ary conditions on our L x L spatial volume is that these constrain us to having an 
odd or even number of vortices respectively. We assume that in both cases we have a 
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Poisson distribution, just as in (l26| ) with A = L"^. Thus the probabihty of vortices 



IS 



where the normahsation factors in the two cases are 



Pin.;L^) = ^r-^-^^ (32) 



E ^^^^^ = sinh(^o^2) (33) 



riv odd 



Mnt= E ^^^^^ = cosh(/.o^2). (34) 



77 ' 

ni, even ^v • 



Let So be the effective action of a single vortex. Then it is a simple calculation to 
obtain 

1 _ {^i,L^Y^ 1 _ {ii,L^Y^ 

= j-^ ^vSo^-^ jrr- E ^^^0 , 

■'^^■^ n„ odd n„ e-uen "'v- 

= So ■ a-^ ln(tanh /i) 

OfJ, 



e 



= 2so ■ . ^— — (35) 

where in the second line we use the shorthand notation yU = ^qL"^ and at the final step 
we have used 2/io = a?a{T) from (|30|). 

In (^) the overall factor of exp{— a^aL^} is no surprise in a linearly confining the- 
ory: it arises from the fact that we are subtracting the partition function for an odd 
number of vortices from that for an even number, just as in (P7|). The prefactor of aL"^ 
arises from the average number of vortices in the volume; both {Stw) and {Snt) are 
proportional to it for large volumes; and hence so is their difference, AS". For small 
enough spatial volumes we expect the periodic volume to have no vortices and the 
twisted volume to have just one vortex. We observe that (|35|) embodies this expecta- 
tion: in the limit — > 0, one obtains AS* — »• sq, since the factor in the denominator 
cancels the cxL^ prefactor in the numerator. Of course, for realistic vortices, the action 
will start to vary once < Ay and this simple formula will break down. 



6 Calculations of AS 

In this section we present the results of our lattice Monte Carlo calculations of AS". 
This is defined (0) as the difference of the average total action on lattices with and 
without a twist: AS = (Stw) — {Snt)- It is related ( [21| ) to the derivative with respect 
to T of the vortex free energy F^. In the low T deconfining phase we expect that 
limi:,^oo -^i; = and hence limL^ooA^" = 0. If we have spatial confinement at high 
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T then we expect this to hold there as welL At low T there are specific expectations 
about how this limit is approached |l|]. At high T we can test the possibilities described 
in the previous section. 

Our individual calculations of {Stw) and {Snt) consist, typically, of 0(10^ — 10'') 
interweaved heat bath and over-relaxed sweeps. We shall begin with a calculation of 
AS in large spatial volumes and for temperatures ranging from well below Tc to well 
above. We shall then present a finite volume study of AS for T > Tc. We shall finish 
with a similar finite volume study of the system in the low T confining phase. 

6.1 AS on large spatial volumes 

Our first calculation is of AS = (Stw) — (Snt) in the limit of large volumes and for 
temperatures ranging from values that are well within the confining phase, to values 
much larger than T^. The vortex free energy (^Tj), and hence AS, should certainly be 
zero for T < Tc, so the interesting question is whether there is a change for T > T^. 
In Fig.|^ we show our results. As is clear, these are compatible with AS = over our 
whole range of (3. However, to assess the significance of this result, there are some 
questions that need to be answered. First, what is the range of T? Then, how large 
is the spatial volume? Finally, and most importantly, are the errors on these values of 
AS small compared to the non-zero value one is attempting to discriminate against? 

First, the range of T. On our Lt = 2 lattice deconfinement occurs for (3 = (3c = 3.47 
(see (p!OD). The range of {3 values for this calculation is /5 G [0.5,9.0]. Using this 
translates into 0.15Tc < T < 2.6Tc. (Note that because of large lattice corrections at 
small /3, the lower limit has to be interpreted with care.) 

Is the volume large? We use L = 12 for /5 < 2.5, L = 16 for 2.5 < (3 < 4.5, and 
L = 24 for (3 > 4.5. First consider T < Tc, i.e. (3 < 3.47. Here a relevant physical 
length scale is provided by ^o- = l/o"(^ = 0). In these units we have aL/^„ > 7.8 
for T < T(.. As T —>■ there is a correlation length, C,p, that diverges (the transition 
is second order) and so there will be some region of (3 near /3c where our volume will 
become too small to accommodate ^p. Such finite volume effects have been studied 
in detail in ^ and one finds that they become significant on a 16^2 lattice only for 
f3 > 3.15 where C,p > 7. Going now to T > Tc we note that there is again an interval 
close to Tc where the correlation length diverges. We do not have direct calculations 
of finite size effects here, but we assume that they are negligible for (3 > 3.7 since 
^p(/3 = 3.7) ~ 3 and ^d{P = 3.7) ~ 6 [Q. In the high-T range our lattice sizes range 
from aL ~ 4:.5^£) at (3 = 4.0 to aL ~ 9^^^ at /? = 9.0. Thus, apart from a window of 
temperatures close to Tc, the volumes we are using are indeed large. (In fact there is 
some reason to believe that even within this window the quantities we are interested in 
will be little affected. For example, the calculated values of the spatial string tension, 
(t(T), show no effects as T passes near Tc. And we see no sign of larger fluctuations 
near Tc in either o"(T) or AS.) We shall later perform some detailed flnite volume 
studies that will conflrm this. 
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How significant is our result that AS is compatible with zero? The extra action of 
a single ultraviolet vortex would arise from the change in sign of a single plaquette 
in each time-slice, and so it would contribute 63^ = Lt{TrUp). This increases from 
SS^ ~ 2.5 at P = 3.0 to 6Sy ~ 3.5 at (3 = 9.0. This is certainly excluded by Fig.|[ 
The action of a more realistic 'fat' vortex will depend on the area, as in (|23|) and (12^) . 
Our results constrain the diameter of such a vortex to be > 10a. Compared to the 
natural length scales in the problem (as discussed above) this is very large, and thus 
we believe that our results for AS strongly discriminate against the twisted system 
differing from the untwisted one by the presence of a single extra vortex. 



6.2 AS as a function of L for T > Tc 

We have seen that limi^oo AS = 0. How is this limit approached? Does it point to a 
gas of confining vortices, and hence a behaviour as in (|35|) ? Or is it just that vortices 
do not exist; for example as in ( pSj) because their size is not stabilised? 

To address this question we have calculated AS* as a function of the spatial lattice 
size, L, at P = 9 and for Lt = 4 which, at this /3, corresponds (|T0|) to T ~ lATc. We 
show our results in Fig.^. We observe that on small volumes AS" does indeed deviate 
substantially from zero. Moreover, we see that the variation is not compatible with 
the AS" oc behaviour expected (p5[) if the 'vortex' spreads to fill the whole spatial 
volume. By contrast it is well fitted by the expression (^) that one derives in a dilute 
vortex gas. (Except for the smallest lattice, which corresponds to a 4^ volume, and 
which can hardly be characterised as a system at a finite T.) 

As an aside we remark that an equally good fit can be obtained using the expression 



in (|35|) without the term in the denominator. Clearly much more accurate calculations 
will be needed to determine whether it is there or not. We have also tried a fit of the 
form (|25|) modified by a screening factor of the form exp(— mL). This gives a rather 
poor fit. 

Returning to the fit using (p5D, we find that the parameter in the exponent is con- 
strained to be 

aV(ar = 0.25) = 0.043 ± 0.002. (36) 

We note that this is in precise agreement with the value of the spatial string tension, 
a'^a{aT = 0.25) = 0.0443(4) that we obtained in (0) from Polyakov loop correlations 
at the same value of a and T. This provides a striking confirmation of the simple vortex 
model. 

Can we learn something about the properties of these vortices? For example, if we 
can estimate the vortex area and also the average density of vortices, then we can 
determine whether the vortices are indeed dilute. The density is straightforward to 
determine. From ( p6| , |30D the average number of vortices per unit area is 

ny = Hq = ^ 0.022 (37) 
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which implies that the typical separation between neighbouring vortices is 



5r = — — ~ 7 (38) 

in lattice units. Now, we have previously suggested that the vortex area should be 
Ay ~ l/{ag^.aT) ~ 9. If so, such a vortex gas will indeed be quite dilute. 

It would be much better to estimate directly; but this is not so easy because of the 
unknown factors, such as Zy{(3). If we ignore these factors then we can obtain wildly 
different estimates of A^. For example, one approach is to take the overall coefficient 
of the fit in Fig.|| 

2soaV = 0.024 ± 0.002 (39) 
and then to use (RBf) to obtain sq ~ 0.27. If we use this value for in (ESf) we obtain 



Ay ~ 4. If instead we use (|29| ) with sq ~ 0.27 then we find A^ ~ 70. This huge 
disagreement suggests that these calculations are much too rough to be useful. At a 
very qualitative level, the fact that (|35|) fits the values in Fig.^ all the way down to 
L = 6, with the same value of the vortex action, does suggest that the vortex is not 
being squeezed in this range. This suggests that the vortex area satisfies Ay < 36, so 
that the vortex gas is at least moderately dilute. 

6.3 AS as a function of L for T < Tc 

In the low T confining phase one expects, from the general analysis of [|l|, that the 
vortex free energy, Fy, should approach zero as L — >■ cxo with a correction that oc 
exp(— a^aL^). This is also the prediction of our simple vortex model if we generalise it 
to a condensate of the arbitrarily long vortex lines that would be needed to disorder 
arbitrarily large Wilson loops at T = 0. 

It would obviously be of interest to extend our high-T calculations of AS* to low T 
so as to test these expectations. In Fig.^ we show the result of such a calculation on 
lattices with = 6 at /3 = 7. Using ([T0|) we see that this corresponds to T ~ 0.7Tc. 
We also show in the figure a best fit of the form in ( P^ to the values for L > 6. This 
clearly fits very well. From the exponent we extract the value 

2 [ 0.043(2) L > 6 
^ " = I 0.048(4) L > 8 

This is nicely consistent with the low-T string tension, a^a ~ 0.0455, that one obtains 
from d^) at /5 = 7. (The comparison with the T = string tension is the correct one to 
make even though T ^ here. The reason is that the spatial string tension does not 
vary significantly with T until T is close to Tc, as we see in Fig.^.) 

As a check of the scaling properties of this behaviour, we have also performed such 
a calculation on = 4 lattices at /5 = 5. Using (^) we see that this corresponds to 
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T ~ 0. 75 Tc. We again show the best fit of the form (pSf ), in this case to L > 4. Again 
this fits very well. From the fits we extract 

, ^ 1 0.091(2) L > 4 

\ 0.080(7) L > 6 ^ ' 

which, we note, is reasonably consistent with the string tension, a^a = 0.0979(13), that 
one obtains by a direct T = calculation [0] at this value of /3. 

So we see that in the low T confining phase, the finite volume corrections to the 
asymptotic AS* = result, are also precisely what one expects from the simple vortex 
model (and indeed from the more general arguments of [Q). We note that the fact 
that the functional form fits all the way down to 4^ and 6^ lattices at /5 = 5 and 7 
respectively, is in contrast to the high-T Lt = 4 case where the value of AS on the 
4^ lattice falls below the fit. The temptation is to interpret this as follows. As L | 
at fixed Lt = 1/aT, the lattice ceases to correspond to a thermodynamic system at a 
fixed T, and so the area of the vortex begins to grow, as it loses its oc 1/T suppression 
factor. This leads to a decrease of Sq and hence of AS*. This is of course no more than 
a plausible speculation. 



7 Conclusions 

That centre vortices might drive confinement is an attractive idea |jT]] which has mo- 
tivated a number of recent lattice calculations, e.g. P, ^. At high temperatures one 
would naively expect to find a gas of such vortices which, in a Euclidean calculation, 
would manifest themselves as vortex world lines that close through the short (= 1/T) 
temporal boundary. In the deconfining phase these should then provide the mechanism 
for 'spatial' linear confinement. The statistical mechanics of such high-T vortices is 
relatively simple, as we saw in a calculation that showed how a spatial string tension, 
(t(T), is generated within a dilute gas of vortices. The observed cr(T) oc T behaviour 
arises naturally if the vortex cross-section has an area oc g^T. 

In contrast to such expectations is the recent suggestion ^ that the vortices might 
not be uncorrelated at large separations, so that there is no linear spatial confinement 
for T > Tc. To address this possibility we performed calculations of the variation with 
T of the vortex free energy and found it to be zero at all temperatures, high and low. 
This supports the conventional expectation of a non-zero spatial string tension. We 
then calculated the approach to zero as the spatial volume increases and found that 
it was controlled by the same spatial string tension, cr{T), that had been obtained in 
earlier direct calculations of that quantity. Moreover the functional behaviour is exactly 
like the one we derived within the dilute vortex gas. From the fits we estimated the 
mean inter-vortex distance and, much more roughly, the vortex area. These estimates 
suggest that the vortices are indeed moderately dilute. 
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We performed similar calculations in the low-T confining phase and showed that 
there too the finite volume corrections are as expected for a dilute gas of vortices, 
although here the argument can be made much more generally 

We carefully reanalysed some earlier calculations H, |^ of the spatial string tension 
and showed that the observed linearity of the potential extends over distances that are 
large compared to the natural dynamical length scales at high T. 

All this provides strong evidence that the usual expectation of an area decay for 
spatial Wilson loops at high T is indeed correct. It also provides support for the 
notion that, at least at high T, vortices do produce this area law behaviour. We intend 
to pursue these questions with a better scaling study as well as a dedicated spatial 
string tension calculation. Such studies are worth pursuing because the dynamics of 
vortices should be much simpler to identify at high T; both in Euclidean calculations 
and possibly through a careful analysis of the relatively tractable D=l + 1 reduced 
theory. 
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Figure 1: The mass of the spatial flux 'tube' versus its length for aT = 1/2 > aT^ 
(□) and for aT = 1/6 < aTc (•). All in lattice units at /3 = 4.5. 
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Figure 2: The spatial string tension as a function of aT at /3 = 9. The deconfining 
temperature, aTc, is indicated. 
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Figure 3: The difference of twisted and untwisted actions on = 2 lattices. T 
increases as /3 increases, and the system deconfines at Pc as indicated. The spatial 
volumes are large in physical units, except very close to P = Pc where there is a 
diverging correlation length. 
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Figure 4: The difference of twisted and untwisted actions at T = l/4a > at /5 = 9, 
for various spatial sizes. Curves correspond to best fits for a dilute vortex gas (solid 
line) and for a single unstable 'vortex' (dashed line). 
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Figure 5: The difference of twisted and untwisted actions at T 
for various spatial sizes. Fit sliown is for a dilute vortex gas. 



l/6a < at /? = 7, 



25 



0.6 



0.4 



^ 0.2 



-0.2 



^ f 



8 10 
L 



12 14 16 



Figure 6: The difference of twisted and untwisted actions at T = l/4a < Tc at /5 = 5, 
for various spatial sizes. Fit sliown is for a dilute vortex gas. 
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